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The frequency-dependent response of a one-dimensional 
fermion system is investigated using Current Density 
Functional Theory (CDFT) within the local approxima- 
tion (LDA). DFT-LDA, and in particular CDFT-LDA, re- 
produces very well the dispersion of the collective exci- 
tations. Unsurprisingly, however, the approximation fails 
for details of the dynamic response for large wavevec- 
tors. 



In particular, we introduce CDFT for the one- 
dimensional spinless fermion model with nearest- 
neighbor interaction, and use CDFT-LDA plus exact 
(Bethe ansatz) results for the groundstate energy as func- 
tion of particle density and boundary phase to determine 
the linear response. The successes and failures of this ap- 
proach are discussed in detail. 



1 Introduction Density Functional Theory (DFT) is 
an efficient and powerful tool for determining the elec- 
tronic structure of solids. While originally developed for 
continuum electron systems with Coulomb interaction IH 
121, DFT has also been applied to lattice models, such as 
the Hubbard model |l3]j4||5J|61, in order to develop new ap- 
proaches to correlated electron systems: lattice models of- 
ten allow for exact solutions which hence can serve as bench- 
marks for assessing the quality of approximations. 

Very useful for applications is the Local Density Ap- 
proximation (LDA) where the exchange-correlation energy 
of the inhomogeneous system under consideration is con- 
structed via a local approximation from the homogeneous 
electron system. A lattice version of LDA has been sug- 
gested for one-dimensional systems Q where the underly- 
ing homogeneous system can be solved using Bethe ansatz. 
For recent applications of Bethe ansatz LDA, see also, for 
example, Refs. ||6l|2l[8l|9l[l0l[II][l2l ■ 

In addition, the time-dependent version of DFT has been 
developed and applied 11311141 . in particular, the current 
density version [|l51 ; a recent review jTSl and book ifTTl 
provide excellent overviews, including the relation to stan- 
dard many-body Green's function approaches. 

In this article, we focus on the one-dimensional spin- 
less fermion model with nearest-neighbor interaction, which 
is exactly solvable in the homogeneous case II19II20I . We 
extend our recent DFT-LDA approach jTS) to current den- 
sity functional theory ifTSl . We present the model in Sect. 
2, and discuss general aspects of linear response theory in 



Sect. 3. CDFT and LDA are presented in Sect. 4; for sim- 
plicity we use the static zero-temperature limit in order to 
present the general ideas. Our results are given in Sect. 5, 
and brief conclusions in Sect. 6. 

In the following, h as well as the lattice constant are 
put equal to one; the system size is denoted by L, and we 
assume periodic boundary conditions. 

2 The model We consider one-dimensional spinless 
fermions described by the Hamiltonian 

H = f + V + Y,vm (1) 

I 

where 

r = -<^(c'^'c+Ci+i+h.c.) (2) 

is the kinetic energy, and 

V^vY,nihi+^ (3) 
I 

the interaction. Generally, the phases {(pi] and local poten- 
tials {vi] can be time-dependent. The hat denotes operator- 
valued quantities. Clearly 

hi = dH/dvi , ]i = dH/d^i , (4) 

where 

i, = -it(c^'^'c+Q+i-h.c.) (5) 
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denotes the local current. Particle conservation follows im- 
mediately by noting that in the Heisenberg picture, we have 



hi = i[H,hi] = -{ji - ji-i) 



(6) 



A gauge transformation is described by the (unitary) oper- 
ator 

U = exp{i^ Xm) (7) 
I 

such that the time development of the transformed wave- 
functions is determined by 



H = UHU+ - xm 



(8) 



Thus a gauge transformation implies the replacements 0; 
4>i~^Xi~Xi+i ™d ^ '^i~Xi- Note that the combinations 
ei = 4>i — {vi+i — vi), corresponding to the electric field in 
electrodynamics, as well as <P — J^i't'i^ corresponding to 
the magnetic flux ||2T1 . are gauge invariant. 

In the context of DFT, we will also introduce an auxil- 
iary single-particle system "s", defined by 



(9) 



where is obtained from T defined in dU by the replace- 
ment (pi (jfi for all I. 

3 Linear response Using standard techniques ll22l 
we discuss briefly the linear response to time-dependent 
perturbations {6vi{ty\, {5(t>i{ty} such that 



5H = ^{hiSvi +ji6(t)i) 



(10) 



The response function x, defined through the relation 



(s{him\ 



/ + 00 / 
dt'x^"\t-t')i 



(11) 

has four entries, Xnn, Xnj, Xjn, and Xjj which each are L x 
L matrices with respect to the lattice sites {Im}; these are 
related to expectation values of commutators, for example 



x'z{t-t') = m-mMt),n„^im 



(12) 



where 0{t — t') is the unit step function, and (...) denotes 
the quantum statistical average. Considering a Fourier trans- 
formation with respect to the time difference, t — t' cu, 
the response functions obey Onsager's relations: 



xZi^;{<Pi}) = xZi^-A-<Pi}) 
xl";{u;;{cf>i})^-xTni^;{^M) 



x';j{io;{M) = x]}\^A-<l>i}) 



(13) 
(14) 
(15) 



The minus sign in ( fT4l i reflects that the current is odd un- 
der time reversal. The density-density, density-current and 
current-current response functions are related to each other 
due to particle conservation. 

Considering a homogeneous situation, i.e., (pi = (f) and 
vi — for all I, we obtain 

X,i„('?,w;(/)) = Xmi(-'?,'^; -(/>) (16) 

Xn]{q,i^;(f>) = -Xjn{-q,i^;-(l>) (17) 
Xjj{'i,^'-,<P) ^ X3]{~q,^^;~(f>) (18) 

where q is the wavevector For the homogeneous single- 
particle system (|9]l the explicit results are as follows: 



Xlf3i<l,^-,(t>) 



nk - nk+q 



L Lu + ek ~ ek+q + iO 



(19) 



where k„„ = 1, = kj^ = and Kjj = w^+^/j- 

In addition, ek = — 2tcos(fc + (p) is the free-particle dis- 
persion, Vk ~ dck I dk the corresponding velocity, and Uk 
denotes the Fermi function. We note also that 

Xnj{q,i^;<P) ^ X37i{q,i^;(f>) (20) 

due to parity symmetry. Particle conservation implies 

■ Xnn = 2 sill(g/2) ■ Xnj (21) 

and 



■ Xnn ^ [2sm{q/2)]^ -Xjj 



(22) 



where we suppressed the arguments {q, ui; 4>) for simplic- 
ity. These relations allow writing the density and current 
response in gauge-invariant form, for example: 



c I . Xnn{q,i^) I . 

2i sm(q/2) 



(23) 



where e{q,uj) = -'iuj(j){q,uj) — 2ism{q/2)v{q,uj); com- 
pare the discussion below (|8]l. Equivalently 



5][q,uj) = e[q,uj) 

such that Lo ■ dn{q, uj) ~ 2 sm{q/2) ■ Sj{q, oj). 



(24) 



4 Current Density Functional Theory In this sec- 
tion, we briefly outline the current density functional the- 
ory and the local density approximation, without discussing 
questions of uniqueness and other mathematical difficul- 
ties ll23l . For simplicity of notation, we restrict ourselves to 
the static zero-temperature limit. The generalization to the 
time-dependent finite-temperature case is straightforward, 
utilizing generating functionals, functional derivatives, etc 
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4.1 General aspects We start with the groundstate 
energy of the Hamiltonian (dJ, E, which is a function of 
the local phases {4>i} and potentials {vi} with the property 



^ (ni) ^ dE/dvi , ji = {ji) ^ dE/dcj)i 



(25) 



We transform to new variables and {ji}, i.e., intro- 
duce the Legendre transform, F, according to the relation 

F = E -Y.i'^^m + ^iji) (26) 



such that 



^dF/dni , (f>i ^ ~dF/dji . (27) 



Obviously, some care will be necessary due to gauge in- 
variance. For the static DFT case, this is a minor problem: 
a constant can be added to the local potentials without es- 
sential changes of the physics. In the general dynamic case, 
one has to keep in mind that density and current are not in- 
dependent, but related by particle conservation. 

In the next step, we perform an analogous Legendre 
transformation for the auxiliary single-particle system (|9]l, 



{vim + cjfiji) 



(28) 



Then, by definition, we have 



Vi = VI + 



dni 



,<Pl^^i + -^T— (29) 



where i^;HXC ^ p _pB^ (jhe superscript "HXC" refers to 
Hartree-Exchange-Correlation.) Introducing explicitly the 
Hartree contribution, E^, through 



^HXC ^ ^ ^XC 



(30) 



with E^ — V J2i "^("^(+1' we arrive at the standard relation 

Vl Vl + Vf + Vr (31) 

where presently up = V{ni+i + the exchange- 

correlation potential is given by vf^ = dE^'~ /dni. In ad- 
dition, 

0^ = 0( + , 0f= = ; (32) 



dji 



here only E^'~^ appears since the Hartree energy does not 
depend on the currents. Explicitly: 



E 



- (o|r+i/|o)-(o,|f^|o,)-£;« + ^0rjz (33) 



where |0) and |0s) are the groundstate wavefunctions of the 
interacting and the single-particle system, respectively. 



4.2 Local approximation In the next step, an approx- 
imation for E^'-^ is needed. As usual, a local approximation 
is employed, according to the following recipe: (i) consider 
the static, homogeneous case, and determine E^'~^{n,j); 
(ii) define e'"'' = E^'^ / L; (iii) approximate jj) 

For step (i), we note that the first term in ( |33] |. hence- 
forth denoted E^^, is known from the Bethe ansatz so- 
lution II19II20I of the model, albeit as a function of n and 
(p. Thus the phase variable has to be eliminated from this 
expression in favor of the current, using the relation j = 
dE^^{n,(j))/d<l> (recall that = L0). In relation to the 
second term in ( |33] |. which we denote E^, we recall the 
single-particle result (— tt/L < <j> < tt/L) 

2t 

E°{(t)) = E°{Q) • COS0 , E^{0) ~ Lsin(7m) (34) 

TT 

where the latter relation holds for large L. Since ~ 
we may expand for small in particular, the Drude weights, 
D^^ and D^, are defined according to the following rela- 
tions (^ 0): 



BA^2 



(35) 

(36) 



Note that D^^ and are functions of the density, and 
D^^ depends on the interaction V. For example, ~ 
(t/7r) sin(n7r) = vp/2tt, where vp is the bare Fermi ve- 
locity, and D^^ = 7rt sin/i/[4/i(7r — /i)] for half filling 
(n = 1/2), where V (in the range ~2t . . .2t) is related to 
fihy V — —2tcosfi. 

Combining the above relations, we obtain 

E^^{n,j) = E^^{n,0) ~ E°{n,0) - LVn^ + ^X'^'f 

(37) 

where 



A-(n) 



1 



1 



(38) 



Note that \'"'{n) < since D^^{n) < D^{n). The next 
steps, (ii) and (iii), are straightforward. The resulting ap- 
proximation may be called CDFT-LDA. 

4.3 CDFT-LDA and linear response For the deter- 
mination of the response functions, we again employ the 
auxiliary single-particle system as follows. First, we con- 
sider the response of the s-system to small variations 6v^, 
Scff, thereby defining the quantity analogous to ( fTTT i. 
(Here and in the following two equations, we will resort 
to a short-hand notation.) Second, we take the variations of 
the Hartree and the exchange contributions to the potentials 
and phases into account, according to 



5v 



pxc 




(39) 
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it is straightforward to verify that 



0.5 



0.2 
0.1 



0.05 



0.02 



27r 



Figure 1 Contour plot of the imaginary part of the dynam- 
ical susceptibility Imx„„(g,ti>) (in units of t^^), as ob- 
tained from (ITTT l. for half-filling and V/t = 1. Here and in 
the following figures, ui is given in units of t. For the plot, 
we replace the imaginary part of the denominator in (fT9] l 
by a finite value, yy, which we choose here to be 0.005t. 



thereby introducing f^xc _ ^ jrxc -pj^g result is 



(40) 



which reduces to the standard RPA expression when f ^'-^ = 
0. Note that ™ only has an nn-entry, in Fourier represen- 
tation given by V{q) = 2Fcos(7. 

The above expression is exact provided the exact func- 
tional f^'-' is used. In the following, however, we rely on 
the results of the previous subsection, and consider in par- 
ticular j 0; in this limit, see above, we may use fj^^*^ — 
0, — 0. Employing particle conservation again, we 

find 



Xrm(g,w) 



fHXC 



4sin^(g/2)^JJ I 



(41) 



where Xnn^ of course, depends on q and lu. If in addition 
tjj^'^ ~ 0, we recover the approximation known as adia- 
batic LDA JTSllBll . 

In order to discuss the result (l4Tl i in more detail, recall 
that in the long-wavelength low-frequency limit the density 
response of the s-system is given by 



(42) 



where the static susceptibility xltat = I/ttv^?. Considering 
the limit lo = Q, q ^ {), and noting that 



fXC 



BA 



2LV 



(43) 



Xstat 



1 a^^BA 

L dn^ 



On the other hand, taking q ^ first, we find 

X,,(g-0,c^^0) = -2i?BA^ 



(44) 



(45) 



i.e., the exact result. The minus sign here is due to our def- 
inition of the response function, compare (fTTT i. Inserting 
(l42T i into iAli . we find (for small q, lo) 



Xnn 



where 



{qv)^ - {uj + i0)2 



1 - ff^'^VF/n Xstat 



(46) 



(47) 



which - as to be expected in view of i45[ - is the exact 
expression. 

The numerical results presented below are based on 
Eqs. dSTll and dTTT i. 

5 Numerical results InFig.[T]we show a contour plot 
of the imaginary part of Xnn {Qi'^) for half-filling and V/t = 
1. The apparent continuum of excitations can be identified 
with the particle-hole continuum; its spectral weight van- 
ishes in the long-wavelength limit. Above the continuum, 
there is a well-defined branch of collective excitations; as 
discussed above, the corresponding velocity for q has 
the exact value. However, the contour plot is almost indis- 
tinguishable from the corresponding one obtained within 
adiabatic LDA (compare Fig. 4 in ifTSl ). 

Thus, in order to highlight the differences between adi- 
abatic LDA ifTSl and the present CDFT-LDA, we present 
in Figs.|2]and[3]lmxTm('?,^) for a fixed wavevector, q = 
7r/2, as a function of frequency, again for half-filling. In 
Fig. |2](^/i = 0.5), even though the interaction is still 
moderate, the spectral weight of the continuum is already 
strongly reduced. However, the frequency range w° . . . uj'^^, 
where uj'^ = 2t\ smq\ and uj^^ = 4:tsm{q/2), is fixed and 
equals the range of the non-interacting model for all V. 
Above the continuum, the collective mode is apparent; for 
this interaction parameter, however, the correction due to 
CDFT-LDA is minimal. 

With increasing interaction, the spectral weight of the 
continuum is further reduced, and the frequency of the col- 
lective mode is shifted slightly to a lower value compared 
to adiabatic LDA, consistent with (l47T i; see Fig. [3] 

The above results have to be contrasted with recent ex- 
act results for the dynamic response of the spinless fermion 
model Il26ll27ll28ll29l . which - unsurprisingly - are not cor- 
rectly reproduced within CDFT-LDA. For example, a con- 
tinuum of collective excitations is found in 02711281 for a 
certain frequency range, < lu < cj^, where uj± are 
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Figure 2 Inix„ri(7r/2, w) versus lo for half-filling, com- 
paring DFT-LDA (green) with CDFT-LDA (blue) for V = 
0.5t. The non-interacting case, V ~ Q (red), is given as 
reference. (77 = O.OOOlt) 







1 




1.5 








e 








1 — 1 








0.5 




— ■ i 


1 i 





1.5 2 2.5 3 3.5 

Figure 3 Same as Fig.lH for V ^ t. 

q- and ;«feracr/o«-dependent. In addition, spectral weight 
is shifted (for a repulsive interaction) to the lower end of 
the continuum, leading to a power-law divergence near uj- . 
For an attractive interaction, on the other hand, a bound 
state is found in the dynamical structure factor above the 
continuum 12811301 . 

In order to clarify how the differences between exact 
and CDFT-LDA results develop with increasing interac- 
tion, we have performed exact diagonalization studies for 
small systems of 16 lattice sites. Some of our results are 
shown in Fig. H] where we plot the imaginary part of the 
local density response function 

Imx"„(t^) = L-'^^lmXnniq,'^) (48) 

for half-filling versus frequency, for V ~ 0, OAt, and 0.8t. 

The usefulness of Im x^ni^) li^^ in the fact that the ex- 
citation energies of the interacting system are given by the 
poles of the response function, which appear as broadened 
(5-peaks in the figures due to the finite value 77 = 0.02t of 
the imaginary part of the frequency. Figure 4a shows the 
susceptibility of the non-interacting system, as reference. 



Obviously exact diagonalization and CDFT-LDA yield iden- 
tical results in this case. 

In Fig. 4b, where V/t ~ 0.4, the two peaks with low- 
est energy are split into doublets. Remarkably, for each 
doublet the position of the peak with the higher energy 
(marked by blue arrows) is almost exactly obtained within 
CDFT-LDA. These peaks correspond to the two lowest pos- 
sible (/-values for a 16 site system, 7r/8 and 7r/4, respec- 
tively, which again demonstrates the validity of the CDFT- 
LDA in the long-wavelength limit. 

On the other hand, the peaks with lower energy within 
the doublets, corresponding to g-values near tt, are clearly 
off. A further feature that is not obtained within CDFT- 
LDA is the appearance of high energy satellites beyond the 
upper limit ui^ = At of the non-interacting continuum. 

The trends already apparent for V/t — 0.4 become 
even clearer for V/t = 0.8 (Fig. 4c). In addition, in the 
exact data more and more spectral weight is shifted down 
to the left sub-peaks of the low energy doublets, a feature 
which is not obtained within CDFT-LDA. The transfer of 
spectral weight to lower frequencies eventually leads to the 
formation of the power-law divergence at the lower end of 
the continuum in the infinite system 1271281 . 

6 Conclusion and outlook We have demonstrated 
that Bethe ansatz current density functional theory correctly 
describes the Luttinger liquid properties of the one-dimen- 
sional spinless fermion model in the long-wavelength low- 
frequency limit, in particular, both limits - u; = Q, q ^ 
and cj ^ 0, (7 = - are recovered; compare (l44l i and ( l47l ). 
The local approximation for the exchange-correlation po- 
tential is insufficient for other aspects: For the static case, 
it misses the "critical" properties related to 2/si?-scattering 
ifTSl . while for the dynamic response, the excitation spec- 
trum for large wavevectors is not correct. It is unclear at 
this moment, at least to us, whether some of the shortcom- 
ings of CDFT-LDA can be cured, for example, by pursuing 
the so-called exact-exchange potential approach lISTI . 

The static case has been discussed in our recent work 
ifTSl in considerable detail. We already noted that the q 
limit of the static response is, by construction, obtained 
exactly within LDA; compare the discussion in connec- 
tion with (l44l) . Furthermore, as a major improvement in 
comparison with the Hartree approximation, Bethe ansatz 
LDA correctly predicts a non-charge-ordered groundstate 
for a large range of parameters. The static density response 
was found to agree very well with the exact result for not 
too large systems, low particle density, and wavevectors 
q < 2kp. 

Good agreement between exact (density matrix renor- 
malization group) calculations and time-dependent DFT is 
also reported in a recent study lfT2l of the collective den- 
sity and spin dynamics of the one-dimensional Hubbard 
model. In this work, the adiabatic local spin-density ap- 
proximation is employed to investigate the density and spin 
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c)V = OM 

Figure 4 Local density response Imx"„(w) versus uj for 
a small (L = 16) system (at half-filling), comparing ex- 
act diagonalization results (green) with CDFT-LDA (red). 
Here 77 = 0.02. The blue arrows indicate excitations that 
are obtained almost exactly within CDFT-LDA, i.e., near 
these peaks the green and the red curves are on top of each 
other on the scale of this plot. 



response to a local time-dependent perturbation for small 
and relatively dilute systems. 

On the other hand, it is obvious that the exact exchange- 
correlation potential can be obtained by numerical meth- 
ods, at least for relatively small model systems. In this con- 



text, we noted previously ifTSl that the Bethe ansatz LDA 
combined with numerically determined exact exchange- 
correlation potentials 1321 might be a useful approach for 
short, inhomogeneous systems, like quantum dots and mo- 
lecules. Ab initio calculations of the linear conductance 
through molecules, e.g., as a function of the gate voltage, 
are usually based on DFT-LDA - but the theoretical and 
experimental results differ typically by an order of magni- 
tude, which seems to be related to the quality of the xc- 
potentials employed; see 1331 for detailed discussions. 

Following this suggestion l32| , we have started an in- 
vestigation of small one-dimensional interacting dots, typ- 
ically consisting of five sites, coupled to one-dimensional 
leads, within (i) LDA, (ii) Hartree-Fock approximation, and 
(iii) DFT with the exchange-correlation potential being de- 
termined from exact diagonalization. Preliminary results 
indicate that the "exact DFT" leads to a considerable im- 
provement compared to LDA and Hartree-Fock; i.e., good 
agreement is obtained with the density matrix renormaliza- 
tion group studies of 1321 . A realistic calculation of trans- 
port properties hence seems to be feasible, by combining 
"exact DFT" for small dots with LDA for the leads, which 
do not necessarily need to be one-dimensional in this ap- 
proach. 
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